Magic ratio of window width to grating period for Van der Waals potential 

measurements using material gratings 
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We report improved precision measurements of the Van der Waals potential strength (C3) for Na 
atoms and a silicon-nitride (SiN^,) surface. We studied diffraction from nano-fabricated gratings with 
a particular "magic" open-fraction that allows us to determine C3 without the need for separate 
measurements of the width of the grating openings. Therefore, finding the magic open-fraction 
improves the precision of C3 measurements. The same effect is demonstrated for a grating with an 
arbitrary open-fraction by rotating it to a particular "magic" angle, yielding C3 — 3.42 ±0.19 eVA 
for Na and a SiNx surface. This precision is sufficient to detect a change in C3 due to a thin metal 
coating on the grating surface. We discuss the contribution to C3 of core electrons and edge effects. 



I. INTRODUCTION 

Van der Waals and Casimir-Polder potentials are the 
dominant interactions between charge-neutral objects at 
nano- to micrometer length scales. As such they have 
attracted considerable interest in the field of microtech- 
nology, as well as the field of quantum-gravity, where 
some theories predict deviations from Newtonian poten- 
tials at very short length scales The dependence of 
these potentials on the geometry and dielectric function 
of surfaces is only partially understood. In the present 
paper we report measurements of sufficient precision to 
study these effects. 

Over the past decade. Van der Waals (VdW) potentials 
between atoms and surfaces have been measured using 
diffraction from nano-fabricated gratings , quantum 
reflection from surfaces [1, and spectroscopy of Cs 
atoms in nano-cells The most accurate of these re- 
ported an uncertainty of 15% @. 

In this paper we report the VdW potential strength 
(C3) for Na atoms and a sihcon-nitride (SiN^;) surface 
with a precision of 5% by studying diffraction of an atom 
beam from a nano-grating. We applied this improved 
precision to show that we can detect a change in C3 due to 
a thin layer of metal deposited on the grating. In addition 
we will discuss the effect of edges on the atom-surface 
potential and we will discuss some of the difficulties in 
calculating the VdW potential for real systems. 

The dominant source of uncertainty in C3 in previ- 
ous diffraction experiments was imprecise knowledge of 
the grating's geometric parameters. Perreault et al. |4|] 
reported an uncertainty in C3 of 25% due to an uncer- 
tainty in the window width {w) of only 1 nm (when 
50 nm). This correlation between C3 and w is the main 
problem that we have overcome by finding the "magic" 
open-fraction. We can now measure C3 independently of 
w. 

Perreault et al. determined w by SEM imaging; 
conventional imaging techniques however cannot easily 
improve these measurements. SEM/ TEM imaging is 
hindered by charging and image charge effects, while 
STM/AFM images show a convolution of the sample and 
the unknown tip shape. The experiment described in this 




FIG. 1: A schematic representation of the grating geome- 
try. The grating bars have a trapezoidal cross-section, with 
a wedge angle a. The grating can be rotated around an axis 
parallel to the grating bars (y) by an angle jS. The period of 
the grating (d) is 100 nm. The width of the grating windows 
(to) varies for different gratings between 40 and 70 nm. The 
thickness of the bars (t) is about 120 nm. w' is the projected 
widow width. 



paper allows us to determine the geometric parameters 
in addition to C3. 

The reason it was difficult to determine either C3 or 
w from a study of diffraction alone is that the effect of a 
small increase in C3 is usually very similar to the effect 
of a small decrease in the window width. In this paper 
we explain the origin of this correlation and we will show 
how the correlation disappears for certain special values 
of the open- fraction w/d. Gratings with this magic open- 
fraction enable us to determine C3 and w independently 
from each other, thus eliminating the effect that limited 
the precision of previous experiments. 

We have experimentally demonstrated this by repeat- 
edly coating a grating with metal until the desired Wmagic 
was reached (Section II). In the study of atom-surface 
interactions however, it is undesirable to contaminate 
the surface. We therefore developed a similar method 
wherein the grating is rotated by an angle, /?,„, to change 
the projected open-fraction such that the C3 - w correla- 
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tion vanishes (Section III). 

Our experimental setup is described elsewhere In 
brief, we used a supersonic beam of Na atoms incident 
on a SiNx material grating with a period of d = 100 nm 
(Fig. [1]). We measured the flux of diffracted atoms as 
a function of position by translating a hot-wire detec- 
tor in the transverse direction (Fig. [5]). We then did a 
least squares fit of these data to determine the diffrac- 
tion order intensities (/„). We measured /„ at different 
grating rotation angles /3, as a function of velocity, rang- 
ing from 1000 to 3000 m/s. In the analysis we used only 
the ratio h/h vs velocity (Fig. ^ to determine both C3 
and w. We choose to focus on these two orders because 
they are more sensitive to C3 than the 0th and 1st orders 
and more easily detectable than higher orders. Studying 
the ratio reduces systematic errors associated with de- 
tector non-linearity and the beam-profile used to fit the 
diffraction data. When measuring diffraction intensities 
as a function of grating rotation, we recorded a complete 
diffraction pattern for each rotation angle, thus taking 
into account the change in separation of the orders due 
to foreshortening of the grating period; this was ignored 
in previous work by Cronin et al. Q. 

Given the typical grating geometry, we know all atoms 
must pass within 30 nm of a surface. We are there- 
fore sensitive exclusively to the short range non-retarded 
VdW potential. We used the model described in refer- 
ence [J| to fit the diffraction intensities. The potential 
was approximated by that of two surfaces of infinite ex- 
tent, CO- located with the two inside surfaces of a grating 
window. The potential of one such a surface is well known 
to have the form 



V{t) 



C3 



(1) 



where r is the distance from the atom to the surface. 
This potential was then assumed to be 'on' only for a 
distance t while the atom is between the grating bars. In 
the WKB approximation, the phase of the atom wave- 
function just beyond the grating is given by 
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FIG. 2: (Color online) Atom flux as a function of detector po- 
sition in the far field, 2.43 m from a rotated diffraction grating. 
Note that at non-normal incidence, the positive orders differ 
in intensity from the negative orders. This asymmetry is well 
described by theory ^ . 
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FIG. 3: (Color online) The ratio of the intensities of the sec- 
ond and third diffraction orders J2//3 is shown as a function 
of velocity for grating Gl at normal incidence. Measurements 
were made after repeated coating with Au/Pd. The curves 
indicate three least squares fits to the data. 
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where z is the direction of propagation and x is the the 
transverse coordinate (see Fig. [1]) . The intensity of the 
far-field diffraction orders can now be calculated using 
techniques from Fourier optics. 



II. MAGIC OPEN-FRACTION 

The question of whether there is an optimal nano- 
grating geometry to use in diffraction experiments has 
been a topic of informal conversation for some time 
We will show that there is indeed a "magic" open-fraction 



that improves the precision of C3 measurements by a fac- 
tor of 5. The VdW-induced phase has curvature, just as 
would be the case for the phase acquired by light due to a 
lens. Thus we see that the VdW potential causes atoms 
to be deflected into higher orders. Similarly, a grating 
with a smaller window width would lead to more inten- 
sity in the higher orders relative to the zeroth order. It is 
tempting to try to approximate the effect of the VdW po- 
tential by a simple diffraction model that uses C^cS = 
and a modified window width 

Cseff = 
WeS = W - Aw(C3). (3) 

This approximation has become the matter of textbooks 
p^ . and it illustrates how a measurement of the VdW 
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FIG. 4: The figure shows contours of = Xmin + 1 in the 
C3-W plane corresponding to the 1-a uncertainty of the pa- 
rameters. Contours are shown for simulated data sets with 
C3 = 3 eVA , d = 100 nm, and different values of w. The 
contours describe extended valleys at most open-fractions but 
narrow to a well defined minimum for one particular magic 
open- fraction. In this example Wmagic = 61 nm, in general this 
value is dependent on C3,q and t, as described in Appendix 
A. 



potential cannot be made unless the window width is well 
known. We derive an expression for Wes in Appendix A. 

The fact that we can make such an approximation 
explains the correlation between the fit-parameters C3 
and w. It is therefore of interest to examine whether or 
not this approximation holds for all C3 and w. We re- 
call that diffraction from a binary amplitude mask pro- 
duces zero intensity in the m-th order when the open- 
fraction w/d = 1/m. However, Cronin et al. Q have 
demonstrated that, in the presence of the VdW poten- 
tial, diffraction from material gratings will never produce 
a diffraction order with zero intensity. We therefore ex- 
pect that there exists a particular magic open-fraction 
where the approximation in eqn. ([S]) breaks down. We 
expect this to occur when the effective open-fraction ap- 
proaches an integer fraction, i.e. 

w Wmagic as Wcs/d 1/m. (4) 

In Appendix A we derive an expression for Wmagic and 
we show that the correlation between C3 and w indeed 
vanishes a.s w ^ Wmagic- 

For simplicity, we first use simulated data to investi- 
gate what happens to the fit-parameters C3 and w when 
w — > Wmagic- Wc simulated I2/I3 vs velocity (similar to 
Fig. 12]) for several grating geometries and fitted the simu- 
lated data with our model. Figure |4] shows a contour-plot 
of the x^-surface in C3 - w space for several different val- 
ues of w. For most values of w these contours describe 
extended valleys, making it impossible to determine ei- 
ther parameter uniquely. Near the magic open-fraction 
however, has a well defined minimum. 

To produce a grating with the magic open-fraction we 
started with a grating that had a larger open-fraction. 



FIG. 5; (Color online) Contours of = Xmin + 1 in C3-10 
space for two gratings (labeled Gl and G2) after repeated 
coatings with Au/Pd. The small contour for "Gl coated 
twice" shows the impact of the magic open-fraction. 



coated it with a thin layer of metal and measured atom 
diffraction. We repeated this procedure until the magic 
open fraction was reached. The samples were coated us- 
ing a Hummer IV sputter coater with a Au/Pd target. 
The coating was applied using a plasma current of 2 mA 
for 30 seconds at a pressure of about 100 mTorr; this 
nominally corresponds to a 1 nm layer of deposition. The 
mean free path length under these sputtering conditions 
is about 1 mm which means that the deposition occurs at 
a wide range of angles (omnidirectional) and that a sim- 
ple geometric model of the atomic trajectories suffices to 
determine how much material is deposited on the inside 
walls of the grating bars. A coating was applied to both 
the front and back of the grating. 

A precise determination of C3 and w requires addi- 
tional knowledge of a grating's geometric parameters. We 
obtain these parameters from rocking curves by doing a 
least squares fit to the data shown in Fig. [HI A careful 
analysis of the covariance matrix for this fit shows that 
there are two parameters, that can be determined inde- 
pendently of C3 and w. These parameters are the wedge 
angle (a), and the rotation angle beyond which no atoms 
can pass through the grating (/3c). 

Armed with this knowledge of a and /3c, we measured 
the ratio I2 / 13 as a function of velocity to determine C3 
and w. These data are shown in Fig. [3l We did this 
several times: after one, two and three coatings for two 
different gratings. Figure [5] shows contour plots cor- 
responding to fits of these data. The shape of the 
contours match very well to the simulation shown in Fig. 
m We can see that for a window width a few nm larger 
or smaller than about 60 nm, the contours describe 

long valleys. The best fit C3 = 4.2 ± 0.4 eVA^ is sig- 
nificantly higher than expected for a SiNj. surface but is 
consistent with a metal coated grating. We will interpret 
the measured values for C3 sections |V] and IVII 
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FIG. 6: (Color online) The intensity of orders 1, 2, 3, and 4 
relative to the 0-th order as a function of the grating rotation 
P (see Fig. [l}. Each set of points at one angle represents 
one diffraction scan hke the one shown in Fig. [2] A least 
squares fit to these data allows us to determine a and /3c 
independently from w and C3. Although w and C3 are left as 
free parameters in this fit, they are not well constrained. The 
data represented here are for a grating determined to have 
a = 5 ± 1 degrees and /3c = 26.5 ± 0.1 degrees. For a further 
discussion of rocking curves see references [1, [O] . 



III. MAGIC ANGLE 

In the study of atom-surface interactions it is usually 
undesirable to cover the surface w\t\v any coatings. We 
can use a technique similar to the one above on a grat- 
ing of any initial open-fraction, by rotating the grating 
around a grating bar. When the grating is rotated by 
an angle greater than the wedge angle a, the projected 
open- fraction is reduced (see Fig. [ij. When we rotate 
the grating, the acquired phase profile changes as the in- 
side surfaces of the grating bar are either rotated towards 
or away from the path of the atoms. Taking into account 
this change, we can rotate the grating by a particular 
"magic" angle /3m such that eqn. (jH) is satisfied. 

A grating was placed on a motorized rotation stage 
with an optical encoder with 1/25 degree precision. We 
acquired diffraction scans at about 50 angles ranging 
from -25 to 25 degrees and fitted the 0-th order inten- 
sity as a function of angle for an a posteriori calibration 
of normal incidence. 

We do not know a priori at what angle to position the 
grating, since we would need to know both the window 
width and C3 to calculate the magic angle. We can how- 
ever detect the magic angle in a diffraction experiment. 
From eqn. (|4]) we expect that the intensity of one of 
the orders will be minimal at the magic angle. Figure 
ini shows the intensities of orders 1 through 4 relative to 
the zeroth order obtained from 50 diffraction scans at 
different angles. The third order is minimal at an angle 
of 10.5 degrees. We therefore expect the magic angle to 
be around 10.5 degrees. 



FIG. 7: (Color online) The ratio I2/ I-i as a function of velocity 
at three different angles of incidence, near the magic angle. 
The solid lines indicate a global fit to the sixteen data points 
that are shown. The data shown here is for the grating labeled 
G3. 



We did this experiment for two different gratings, one 
that was pure SiN^; and one that had been coated with 
Au/Pd. The results for the Au/Pd coated grating are 
shown in Fig. [7] and [S] The magic angle was found to 
occur at 10.5 degrees. To achieve the best precision in C3 
and w we did a global fit to hlh vs velocity at 9, 10 and 
11 degrees, shown in Fig. [71 The combination of C3 and 
w obtained this way is consistent with diffraction data at 
normal incidence, as can be seen from Fig. [Sj We find 

C3 = 4.8 ± 0.5 eVA , which is once again significantly 
larger than the value expected for a SiNa; surface. 

The same experiment on a grating with no Au/Pd 

o 3 

coating on the other hand yields C3 = 3.26 ± 0.16 eVA , 
which is consistent with theory for a SiNj; surface. The 
contour for the clean grating is shown in Fig. [HI The 
values found for C3 for different gratings and coatings 
are summarized in Table U] The fit parameters for both 
gratings are summarized in Table HIl In order to find the 
physical C3 for this combination of atom and surface, we 
need to take into account the shape of the potential more 
carefully. This is discussed in section IVl 

The grating labeled G3 was used by Lepoutre et al. 
[1^ to measure VdW induced phase shifts in an atom 
interferometer. Our determination of the grating geom- 
etry enabled Lepoutre et al. to verify the form of 
the VdW potential. 

IV. THEORETICAL PREDICTIONS 

An ab-initio calculation of the the atom-surface inter- 
action can be a very challenging undertaking; even a cal- 
culation for idealized atoms and surfaces can only be done 
analytically in a hand-full of special cases. For an arbi- 
trary geometry there is as of yet no proven method to 
exactly calculate the interaction strength [l3l| . For a real 
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and the surface of an ideal conductor as compared to a 
model not including core electrons. For a SiN^; or a gold 
surface, theory predicts a more modest increase in C3 on 
the order of 6%. This is understandable since core exci- 
tations mainly contribute to apoi(ii-i-') at high frequencies 
where the response of a real surface is small. Our mea- 
surements using a SiN^; surface are in better agreement 
with the model that includes core-electrons, though the 
deviation from a simple Lorenz oscillator model is only 1 
standard deviation. The results are summarized in Table 
III 

V. PAIRWISE INTERACTION 



FIG. 8; (Color online) Contours of — Xmin 4- 1 in C3-W 
space for the Au/Pd coated grating labeled G3. The blue 
(thin solid) ellipse corresponds to a fit of a single diffraction 
scan like Fig. [5] at 2000 m/s at normal incidence. The black 
(dashed) and green (dotted) ellipses correspond to data sim- 
ilar to Fig. [7] at a single angle /3 = and P = 10.5 degrees 
respectively. The red (bold solid) ellipse corresponds to the 
global fit of the three angles shown in Fig. [7] resulting in the 
best constraint on C3. 



system, we must consider the frequency response of the 
atom as well as that of the surface, which is affected by 
the composition and the geometry of the surface. For 
an infinite plane surface, the VdW coefRcient in the non- 
retarded regime is given in reference |T3 [ by 



C3 



h 

47r 



e zo; - 1 
apoi(«w — — — — - duj 5) 
e ^'^j + 1 



where apoiii^) is the atomic polarizability and e{iuj) is 
the dielectric response function of the surface in atomic 
units. In general apoi(ii^) is given by 



apoi(«w) 



fn 



(6) 



where /„ are the oscillator strengths for transitions from 
the ground state to all the n other states with ^ fi = 1. 
For the present case it is sufficient to include only the two 
D-lines of sodium since their combined oscillator strength 
foi + /d2 = 0.961. Indeed, for w = eqn. ^ yields 

Q!poi(0) = 160.7 a.u. = 23.81 which accounts for 99% 
of the static dipole polarizability of sodium. 

The optical response of the silicon nitride was obtained 
experimentally in reference P^]. The material studied 
in this reference was produced in a way similar to the 
material of our gratings. The stoichiometry of the SiN^; 
is not exact, but x is nearly equal to 4/3. 

Theoretical calculations by Derevianko et al. [l^ have 
shown that the effect of core electrons to C3 can be sig- 
nificant. They report a 15% increase in C3 for sodium 



The parameter C3 is only defined for an atom near 
an infinitely extended plane surface for which the poten- 
tial has a well known C^/r^ dependence. In the case of 
a more complicated geometry such as our gratings, the 
potential is in principle a function of all three spatial 
coordinates (VVdw = Vydwix,y, z)). There are multiple 
ways of defining a C3 in this case. One way is the ap- 
proximation discussed before eqn. ([T]), where we use the 
potential of two infinite planes co-located with the grat- 
ing's inner walls. This approximation however ignores 
the effect of edges. A more precise way to define C3 is to 
use the pairwise interaction (PWI) model. 

The PWI approximation assumes that the interaction 
of an atom with a solid body is proportional to the sum 
of the interactions with each of the atoms composing the 
body. Mostepanenko et al. suggested that the PWI 
model, which neglects any multi-body interactions and 
screening effects, gives the right spatial dependence of 
the potential, but requires a normalization constant in 
order to yield the correct potential [13. The PWI po- 
tential is given by 



V{r) = [ dV'- 
Jv 



nCe 



„/|6 



(7) 



where the normalization constant K is dependent on the 
geometry of the solid, and its material properties, n is 
the number-density of atoms and Cg is the atom-atom 
interaction constant. 

The normalization constant can be obtained by com- 
paring the PWI model with an exact calculation for an 
infinite plane. We can then express KnC'e in terms of C3 
for an infinite plane: 



V = - dx dy 

Jx>Q J 
_ -f^planoTTnCg 



dz 



plane 
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where d is the distance to the surface. We can now iden- 
tify C3 = TrKpis^ncnCe/d. A numerical evaluation of the 
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integral in eqn. ([2]) using the potential in eqn. ([T]) with 
KuCq — GCs/tt, shows that the the PWI model yields a 
smaller 4'Vdw by about 5% as compared to the approxi- 
mation discussed before eqn. We accordingly apply 
a 5% correction to the C3 found from our fits. 

By construction the PWI model gives the correct po- 
tential near the middle of a grating bar and close to the 
surface, where the effect of edges is negligible. Near an 
edge, the PWI model may yield a potential that differs 
from the actual potential by no more than 13%. We ob- 
tained this upper bound to the error by comparing the 
PWI model to the exact solution in the limiting case 
of an atom far away from a perfectly conducting sphere 
fisl [T9I I . In this case the PWI potential is expected to dif- 
fer most from the real potential by reasoning analogous 
to that in reference 17] . We further assume that the real 
potential of an atom near a surface of finite extent must 
be smaller than the potential near an infinite surface thus 
creating an additional upper bound to the error. Based 
on an analysis of the error along the propagation path 
of the atom, we conclude that the PWI correction intro- 
duces an error that is no larger than 3% to the acquired 
phase, and thus to our measurement of C3. 

VI. THE EFFECT OF A THIN LAYER OF 
METAL 

The two measurements of C3 for Na and an Au/Pd 
coated grating both showed a significantly larger C3 than 
is expected for a clean SiN^, surface. The magic open- 
fraction experiment and the magic angle experiment were 

o 3 

done using two different samples and give C3 = 4.3 eVA 

o 3 

and C3 = 4.8 eVA respectively. These results are in 
agreement with theoretical predictions for C3 near a gold 
coated surface. 

A thin surface layer (roughly less than an optical skin- 
depth) does not produce the same VdW potential as a 
bulk material, even in the idealized case of a uniform layer 
that retains the optical response properties of the bulk 
material. Studies of the Casimir potential between large 
bodies have detected a dependence on surface properties 
at distances greater than 100 nm |23|. For thin layers, less 
than a skin-depth of the surface material, the effect of the 
surface layer is significantly reduced at these distances. 
Based on theoretical work however, we expect the effect 
on the short range VdW potential to be significant even 
for very thin layers [22] . 

Our metal coated gratings consist of a thin coating of 
a Au/Pd mixture on top of a silicon nitride nano-bar. In 
modeling such a system we must take into account the 
interface between the vacuum and the thin outer surface 
as well as the interface between the surface layer and the 
underlying bulk material. To do this, we have evaluated 
expression 4.14 in reference [22j. We used a Drude model 
for eAu(«w) and an insulator model for esiN^(«w). 

Two problems present themselves in a precise interpre- 
tation of the results. First, the atom-surface potential is 



dependent on the thickness of the coating, which is not 
very well known and may be non-uniform. Second, in 
the presence of a thin surface layer the VdW potential 
no longer follows an exact potential. For a 2 nm 

□ 3 o 3 

layer of Au, C3 varies from 4.8 eVA to 4.2 eVA at dis- 
tances of 2 nm to 25 nm from the surface respectively. 
Where we used C3 = V{r) -r^ We evaluated the potential 
at a distance of 10 nm, where our experiments are most 

sensitive A 2 nm layer of Au on the surface gives 

= 3 

C3 = 4.5 eVA . This value depends on the thickness of 
the Au layer, a 3 nm layer for example gives C3 = 4.6 

eVA'. 

In the analysis of our data, we still assumed a 1/r^ 
potential since we are most sensitive to a small range of 
distances and the physical potential only deviates weakly 
from this form over a short range. Any errors introduced 
by this assumption must be smaller than the range of 
C3 given above and are of the order of the uncertainty 
reported. 

Despite these difficulties, the reported C3 for a metal 
coated grating is significantly different from the SiNx 
measurement. The change of more than 30% is consis- 
tent with theory (see Table Our measurements on 
both metal coated samples are clearly inconsistent with 
a clean SiNa; surface. 

The growth of thin films of metal is a complicated pro- 
cess subject to many experimental parameters [IJ] , how- 
ever there are several pieces of evidence that show that 
there is a significant amount of metal on the inside sur- 
faces of the grating. First, we know that the mean free 
path length in the coating process is about 1 mm. This 
means the deposition is omnidirectional and the inside 
of the bars will be coated with nominally 1 nm of metal. 
Second, SEM images of the metal coated grating show 
significantly reduced effects of charging as compared to 
before coating, this indicates the surface is conducting 
charge. Third, x-ray spectra we obtained in the SEM 
show the presence of Au. Finally, the measured C3 is 
consistent with metal coating on the inside of the bars. It 
is however inconsistent with both theoretical predictions 
and our measurement of C3 for a clean SiNj, grating. 



VII. THE EFFECT OF THE NA BEAM ON THE 
GRATING 

During an experiment, a grating is exposed to the atom 
beam for several hours. A set of 50 diffraction scans used 
to make Fig. [S]takes about 30 minutes to acquire. Such a 
data-set is taken at four different velocities to fully char- 
acterize a grating. The total exposure time to the atom 
beam is about 2 hours per experiment. To determine 
if prolonged beam exposure causes contamination and a 
affects C3, we exposed the grating to a high flux atom 
beam. Under normal experimental conditions the atom 
beam is coUimated by two 10 fj-m wide slits which atten- 
uate the flux at the location of the grating by a factor 
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TABLE I: Experimental results and theoretical predictions 



Experiment 


Ca (eVA') 


This work: 






Na and SiN^ fit parameter, grating A 


3.26 


± 0.16 " 


Na and Au/Pd fit parameter, grating Gl 


4.30 


± 0.5 


Na and Au/Pd fit parameter, grating G3 


4.80 


± 0.5 


Na and SiN^ with PWI correction, grating A 


3.42 


± 0.19 " 


Na and Au/Pd with PWI correction, grating Gl 


4.51 


± 0.5 


Na and Au/Pd with PWI correction, grating G3 


5.04 


± 0.5 


Previous work: 






Na and SiN^ [4] 


2.70 


± 0.8 


Theory 




Na and perfect conductor 


7.6 


[16] 


Na and SiNa, using single oscillator Qpoi(ia;) and tabulated e(ia;) 


3.3 


[15], eqn. (O 


Na and SiNj; using tabulated apoi(ii:i;) and e{iijj) 


3.48 


[15, 2^, eqn. © 


Na and Bulk Au 


5.11 




Na and 1 nm Au at r = 10 nm 


4.3 


[22] 


Na and 2 nm Au at r = 10 nm 


4.5 


[22] 


Na and 3 nm Au at r = 10 nm 


4.6 


[22] 



"The sources of the reported uncertainty are discussed in Ap- 
pendix B. 



of about 100. To test the effect of beam exposure we 
removed both coUimating slits and exposed the grating 
to the unattenuated atom beam source for a duration of 
2 hours. 

The results are shown in Fig. [9] After the first 2 
hours, the open-fraction is not significantly affected, after 
another 2 hour exposure, the window width is changed by 
2.5 nm. Table HIl shows the best fit geometric parameters 
before and after coating, Aq represents a clean grating, 
Ai and A2 correspond to the same grating after one and 
two exposures respectively. 

The best fit value for C3 is not significantly changed 
after exposure to the atom beam. We propose this is due 
to non-uniform coating of the grating. The directional- 
ity of the atom beam suggest that the majority of the 
atoms would hit the front face of the grating bars, only 
about 10% would hit the side of the bars directly. AFM 
images confirm that the sodium forms large clumps on 
the front face of the grating that overshadow the grat- 
ing windows. The change in the wedge angle after the 
coating is consistent with this explanation. 

We conclude that under normal experimental condi- 
tions, the effect of beam exposure should be 100 times 
smaller corresponding to 0.01 nm per 2-hour experiment. 
These deposition rates are consistent with the atom count 
rate we measure at the detector. A small amount of con- 
tamination will still be present. 




After 1x2 Hour Exposure 
Clean Grating 



After 2x2 Hour Exposure 



2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 
CjCeVA') 

FIG. 9: The effect of exposure to the Na atom beam is demon- 
strated. One 2 hour exposure in this figure corresponds to a 
dose 100 times larger than in a typical experimental run. The 
clean SiN^ grating shows a C3 of 3.26 ± 0.16 eVA 



VIII. CONCLUSION 

We measured the VdW potential for Na and a SiN^; 
surface with 5 times better precision than previous 
work. We made use of the fact that using a grating with 
a particular magic open-fraction or a magic rotation 
angle, C3 can be determined independently from the 
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TABLE II: Best fit parameters found from least squares fits 
for three different gratings. The origin of the quoted uncer- 
tainties is explained in Appendix B. Ao represents a clean 
grating, Ai and A2 correspond to the same grating after one 
and two exposures to the Na beam respectively. Gl and G3 
are Au/Pd coated gratings. G3 is the Au/Pd coated interac- 
tion grating used in reference [l^ . 



with 



Grating 


C3 


(eVA^) 


w 


(nm) 




(deg) 


t 


(nm) 


name 


















Ao'"= 


3.26 


± 0.16 


57.0 


± 1.0 


3.5 


± 0.5 


140 


± 10 




3.24 


± 0.3 


57.5 


± 1.0 


3.3 


± 0.5 


140 


± 10 


A2'"' 


3.1 


± 0.3 


54.5 


± 1.0 


3.1 


± 0.5 


140 


± 10 


Gl''= 


4.3 


± 0.5 


61.6 


± 1.0 


4.5 


± 0.5 


110 


± 10 


G3"'= 


4.8 


± 0.5 


53.0 


± 1.2 


5.0 


± 1.0 


110 


± 10 



"Ca and w determined by magic angle method. 
'C3 and w determined by magic open-fraction method. 
■^Clean SiN^, grating. 
■^Grating exposed to Na atom beam. 
Au/Pd coated grating 



geometric grating parameters. This method also yields 
a precise determination of the geometric parameters 
of the grating. These measurements are not subject 
to the systematic problems that plague conventional 
imaging techniques and are therefore useful in other 
experiments using nano-gratings. Our measurements are 
precise enough to detect an increase in the atom-surface 
potential due to a thin layer of metal. We also detected 
the effect of extended exposure to the atom beam on the 
grating parameters, however the effect of contamination 
by the atom beam was verified to be negligible under 
normal experimental conditions. Our measurements 
are now approaching a precision where they are sensi- 
tive to edge-effects and the contribution of core electrons. 

This material is based upon work supported by 
the National Science Foundation under Grant No. 
PHY-0653623 and the Arizona TRIP imaging program. 
The authors thank the members of the Vigue team in 
Toulouse and Oleg Kornilov and J. Peter Toennies for 
many fruitful discussions. 



Appendix A: Deriving the magic open-fraction 



C3t 

hv 



{w/2~xf {w/2 + xf 



(A2) 



where w^d and t are the grating window width, period 
and thickness, v is the velocity of the atoms and x is the 
position along the grating window (Pig. [Ij. This integral 
is represented graphically in Pig. [lU] by a Cornu spiral 
0, [2^ ■ The length of the curve is equal to the window 
width, while the distance between the endpoints gives 
Por C3 = the curve would lie on the circle shown 
in the figure; when C3 is non-zero the extremities of this 
curve spiral away from the circle, one inside the circle 
and one outside. This immediately demonstrates that in 
the presence of the VdW potential, there are no orders 
with zero intensity. 

To study the effect of small changes AC3 and Aw on 
the diffraction intensity, we need to look at the corre- 
sponding translation of the endpoints. As indicated in 
the figure, a change in w adds a length in the middle 
of the curve where the curvature of the Cornu spiral is 
nearly equal to that of the circle, thus causing the end- 
points of the spiral to move parallel to the circle. A 
change in C3 increases the curvature of the end of the 
spiral, which has an effect on the endpoints that is simi- 
lar but not exactly the same. 

To equate a change in C3 to an equivalent change in w, 
we consider only the projection of the endpoint onto the 
circle that corresponds to C3 = 0. This neglects radial 
translations of the endpoints due to AC3. The points 
on this circle nearest to the end points we label ±Wcff/2 
for reasons that will soon become clear. To approximate 
the location of these two points we can use the fact that 
the part of the spiral where (j^vdw > can be neglected. 
At the point on the Cornu spiral where (fivdw = ""j the 
tangent line of the spiral is antiparallel to the circle. The 
value of WcB can then be found from eqn. (jAip by solving 



exp 



2?™ , 

in — i X 

a 



-exp 



27m 
-i——Wcs 2 
a 



(A3) 



where x' is the coordinate where (f>vdw — tt. The minus 
sign on the right hand side accounts for the fact that the 
tangent vectors on the circle (C3 = 0) and on the Cornu 
spiral are pointed in opposite directions. The value of 
x' can easily be found by considering the potential of a 
single wall, it is however more convenient to define xq ~ 
w/2 — x' , the distance from the atom to the nearest wall. 
We find xo by solving 



To better understand the correlation between the win- 
dow width and the VdW coefficient, we consider diffrac- 
tion from a grating with rectangular bars. In the far field, 
the intensity of the n-th order is /„ = Where the 

complex amplitude An is shown in reference [4| to be 



±9l 

TlV Xq 



(A4) 



The relationship between Aw and AC3 is then seen to 
be 



1 



0/2 



^ri = -7 / exp 

" J-w/2 



2'Knx 



+ ipvdwix) 



dx (Al) 



Aw = — — AC3 
0C3 



(A5) 
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0.20- 



0.15- 



- 0.10- 




-0.05 0.00 0.05 0.10 

Im ( S ) 

FIG. 10: (Color online) A Cornu spiral is a graphical rep- 
resentation of the integral in eqn. (|A1|I . it is a curve in 
the complex plane parameterized by {^[^(g)], where 
Sil) = lfL/2^''P[H~^+^vdw{x))] dx. The length 
of the curve is equal to the domain of integration (w), while 
the distance between the endpoints gives the absolute value 
of the amplitude of the diffraction order. The figure shows a 
Cornu spiral in blue (dark grey) for the second order (n=2) 

° 3 

for a grating with w = 50 nm and C3 = 3 eVA . The fig- 
ure also shows a Cornu spiral for a grating with the magic 
open- fraction in red (light grey). 



Using common values for the parameters, t = 120 nm, 
V = 1000 m/s and C3 = 3 eVA^ we get 



AC3 
Aw 



nhv 



-1/3 



0.8 



eVA' 



(A6) 



This is very close to the relationship found empirically 
in reference [3] between the best fit value for C3 and the 
fixed value used for w. 

A rough approximation for /„ can now be found by set- 
ting C3 — and using Wcs in place of the physical value 
for the window width. This approximation is not very 
useful in interpreting experimental data but it does allow 
us to make an estimate of the magic open-fraction. This 
approximation predicts a zero in the single slit diffrac- 
tion envelope at order n when w/d = m/n for integer m. 
When {n, m} = {2,1} the second order is predicted to be 
missing. Since there are no diffraction orders with zero 
intensity, we find the magic open-fraction by solving 



and C3 indeed become small near w = 60 nm. There are 
also magic open-fractions near w = 45 nm and w — 90nm 
corresponding to {n, m} = {3, 1} and {2, 2}. In general 



^magic 



n \TrTiv 



1/3 



(A8) 



Figure [TU] also shows a spiral for a grating with the 
magic open-fraction, demonstrating how the model using 
only WcS gives the wrong result. 

Near the magic open fraction, the approximation of 
using an effective open fraction with C3 = is no longer 
valid. There is no longer a simple analytic expression for 
the relationship between the fit parameters C3 and w. 
In order to still quantitatively explore AC3/ Aw near the 
magic open fraction, we have numerically computed this 
derivative. 

We can now show that the correlation between the fit 
parameters C3 and w indeed vanishes (i.e. AC3/AW = 0) 
at this magic open-fraction. For a single measurement 
at a single velocity /„(ui), we find a relationship be- 
tween C3 and w that is, at least locally, linear with slope 
ai = dcsln/dwin- This relationship is shown in Figure 
[TTlbv the solid line. If we include the uncertainty associ- 
ated with /„ the linear relationship becomes a confidence 
region, given by the area between the dashed lines. In 
order to determine the combination {C3, w} we need an- 
other measurement at a different velocity /„(u2). If the 
slope ai is similar to 02 the confidence region is still an 
extended valley. Only if the angle between the two lines 
is large, are both parameters C3 and w well constrained. 
In that C&SG, clS indicated in Fig. [TTJ X marks the spot. 
From Fig. [TT]we can also see that the condition of mini- 
mum correlation is not exactly the same as the condition 
for the best constraint on w and C3. The former only 
requires ai = —02 while for the latter we want the an- 
gle between the two lines to be 90 degrees [26|. The 
fact that dcsln/dwin can be both positive or negative 
means that a suitable combination of measurements can 
yield AC3/AW = Cov(C3, w)/Var(u;) = for the two 
fit-parameters. This quantity can be also be obtained by 
fitting a set of data with w fixed and determining C3 for 
different values of w [13] . 

The angle between two lines with slopes ai and 02 is 
given by 7 = atan(ai) — atan(a2). Figure fT2l shows 7 as a 
function of w and shows that the angle indeed approaches 
90 degrees around the predicted open-fraction Wmagic/t^- 
The largest 7 occurs when one of the slopes is negative. 
This only happens in a narrow range near Wmagic, when 
either dc-^In or dwin changes sign. We therefore expect 
the best constraint on C3 to be obtained with an open- 
fraction (or grating rotation angle) that minimizes /„ , as 
our simple guess predicted. 



Woff/ d = (i 



'magic 



2xo)/d^ 1/2 



(A7) 



Using C3 w 3 eVA , we find 2xa ~ 10 nm so we predict a 
magic open-fraction for w = nm. A numerical calcu- 
lation shows that the covariances of the fit parameters w 



Appendix B: Discussion of reported uncertainties 

Given the small uncertainty in the reported values for 
C3, it is appropriate to discuss the various sources of 
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60.0 ■ 
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2.6 



— I — 
2.8 



— I — 
3.2 



— I — 
3.4 



— I — 

3.6 



FIG. 11: (Color online) The two solid lines correspond to the 
relationship between the best fit C3 and w for a measure- 
ment of I2 at velocities vi = 1000 m/s and V2 =2000 m/s re- 
spectively. The dashed lines correspond to the la confidence 
interval. The area where the confidence intervals overlap is 
smallest when the angle 7 between the two lines is closest to 
90 degrees. 




58 60 62 64 66 68 

w (nm) 

FIG. 12: The angle 7 between the two lines in Figure [11] is 
plotted versus w. The angle is closest to tt/2 around w — 61 
nm which is indeed around where we expect the magic open- 
fraction to be. As the angle approaches n/2 the parameters 
C3 and w become maximally constrained. 7 = n/2 is indi- 
cated by the dashed line. 



statistical and systematic errors. Errors due to surface 
roughness or variations in the grating parameters have 
previously been discussed in the literature. First, we con- 



sider variations in the grating parameters such as the bar 
thickness. SEM images that show that the RMS devia- 
tion of the grating bar width is about 1.4 nm. We do not 
observe the exponential dampening predicted in reference 
[2^ so we think surface variations can not be described 
in terms of a distribution of the geometric parameters. 
This is consistent with findings in reference [3|. Instead 
they can be treated as surface roughness. 

In quantum reflection experiments in particular, sur- 
face roughness influences the reflection efficiency. In our 
experiment however, surface variations on a scale smaller 
than the grating bar thickness tend to average out since 
the atoms sample multiple areas of the surface. The po- 
tential near a rough surface, to first order has the same 
spatial dependence and the same C3 as a smooth surface 
but relative to an effective surface distance [29|]. This 
should be taken into account in the interpretation of our 
report of geometric parameters, but is of less significance 
to our reported values of C3. 

We determine the velocity of our atom beam by mea- 
suring the distance between diffraction orders. The un- 
certainty in the velocity is due to the uncertainty in the 
distance between the grating and the detector which we 
know to within 1%. 

As can be seen from Table IIIIl the dominant source of 
uncertainty is the uncertainty in determining the 12/13- 
This is due to the fact that one of the orders is very small 
near the magic open-fraction or magic angle. The total 
error given by the individual errors added in quadrature, 
is 5.7%. 

The PWI correction factor of 5% was obtained in a 
simplified geometry (using bars with rectangular cross- 
section) by computing the correction to the integrated 
phase (eqn. ^) using the PWI potential. We did this for 
multiple paths through the grating at normal incidence. 
Since this is a computationally intensive procedure it was 
not feasible to fit our experimental data with this model 
and get more exact knowledge of the effect of the PWI 
correction on the fit-parameters. 

The main source of uncertainty in the determination of 
w is the uncertainty in the wedge angle a. The contours 
shown in Figures [5] and [S] show a far more stringent bound 
on w than we ultimately quote in Table |lT] because the 
X^-contours correspond to a fit with a fixed. We note 
that the parameter w refers to the average window width 
so even though we allow for variations in the physical 
window width from point to point on the order of 2 nm, 
the parameter w can be constrained to within 1 nm. 
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TABLE III: Sources of errors for measurement of C3 for a Na atom and a SiN^, surface. 



Quantity 


Value 


Uncertainty 


Uncertainty in C3 


Grating-detector distance" 


2.43 m 


1% 


1.0% 


Grating wedge angle a 


3.5 deg 


0.5 deg 


0.9% 


Critical angle Pc 


26.5 deg 


0.1 deg 


1.5% 


Geometry variations 


2 nm RMS max 


2 nm 


1.0% 


Surface roughness 


2 nm RMS max 


2 nm 


1.0% 


I2/IA from fit of diffraction data 


1 - 7 


10% 


4.2% 


Total uncertainty in fit parameter C3 


3.26 eVA'" 


0.16 eVA'' 


4.9% 


PWI correction 


5% shift 


3% 


3.0% 


Physical C3 
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"'The grating-detector distance is required to determine the veloc- 
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